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Models of four delta wings were built and tested in the 8 
foot by 8 foot transonic wind tunnel at the National Transonic 
Facility, NASA Langley Research Center. The wings are identical 
in planform shape with a swept-back angle of 65 degrees, but bear 
different leading edge profiles. The models were tested under 
pressurized and cryogenic conditions to simulate true flight 
Reynold's numbers. Data on the aerodyanamic forces and pressure 
distributions at various locations on the wings were taken at 
various flight Mach numbers and angles of attack. Effects of high 
Reynold's number and leading edge radius on the aerodynamic 
characteristics of the wings are being accessed. 

A preliminary analytical study was done at MIT, but it was 
limited to the case of low Reynold's number, thin layer laminar 
flow over one of the wings [1]. To thoroughly understand the 
turbulent, vortical flow around the wings, an effort to perform 
computational aerodynamic analysis is being committed. The 
objective of the analysis is to supplement and validate the 
experimental data and explain the high Reynold's number and 
leading edge effects. GRIDGEN, a software developed by General 
Dynamics, is being used to generate the grid topology for the flow 
field around the wings. The flow solver to be used is CFL3D, a 
computational fluid dynamics code developed by the Computational 
Aerodynamics Branch at Langley. 

Based on the geometric description of the wings [2] , a 
Fortran program called WINGSURF has been written to generate the 
databases defining the surface geometry of the wings. To match 
the true geometry of the models in the wind tunnel for realistic 
comparision with experimental results, databases defining the 
sting support for the wing models have also been created by two 
other Fortran programs, STINJOIN and STINREAR. Listing of the 
programs are attached and the geometry of a typical wing model 
with the sting support is shown in Fig. 1. 

From the geometric databases, the computational grid for each 
wing both with and without the sting was created interactively on 
an IRIS Graphics Workstation employing the single block approach. 

A computational domain of 20 times the chord length of the wing in 
all directions was taken. Figs. 2 and 3 show the grid topology on 
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* wing surface and in the entire computational domain 

respectively, for free flight (no sting) condition. Figs 4 and 5 
V topology for a typical wing model attlch„? ?S 
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p ints. Also, changing the leading edge profile would reouire a 
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computational domain. This means excessive computing time. 
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computational domain is divided into seven blocks, as shown in 
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Surface grids for each block are in the process of be'inc 
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computational domain. The flow solver can then be applied to 
compute flow properties in the domain. Results wil^be used to 
compare with wind tunnel data and to access effects of hich 
Reynold s number and leading edge radius. g 
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program wingsurf 

* dimension a(4),b(4),c(4),d(4),x(22,22),z(22,22),y(22,22) 

pi = 4.*atan(l. ) 
del = 65./180.*pi 
„ rc = 12.*tan(del) 

xle = . 15*rc*cos(del) 
xte = . l*rc 

data bl/. 510024/, cl/-. 19819008/, dl/. 025671542/ 

" data a/. 4935262, .34897572, .20148123, 0.0/ 

data b/. 080048085, .29224873, .5087712, .804546399243/ 

data c/-. 19710469, -.28382228, -.37230603, .49317707317/ 

- data d/. 046323876, .062270923, .078542758, .100770498643/ 

do 900 k=l,4 
do 100 j=2,22 
xj = (22-j)/20. 

xn = xle*(exp(5.*xj)-l. )/(exp(5. )-l. ) 
x(l,j) = xn/cos(del) 

yd»J) = 0. 

z(l,j) = a(k)*xn**.5 + b(k)*xn + c(k)*xn**2 + d(k)*xn**3 
100 continue 

z(l,2) = bl*xte + cl*xte**2 + dl*xte**3 

x(l,l) = . 15*rc 

x(2, 1) = . 2*rc 

x(3, 1) = .4*rc 

x(4,l) = .6*rc 

x(5, 1) = rc - 9.395 

x(6,l) = . 8*rc 

x(7,l) = . 9*rc 

do 110 i=l,7 

y(i,l) = 0. 

z(i,l) = z(l,2) 

110 continue 

do 130 i=8 , 22 
xi = (22-D/15. 

xi = xte*(exp(4.*xi)-l. )/(exp(4. )-l. ) 
x(i,l) = rc - xi 

y(i,l) = 0. 

z(i, 1) = bl*xi + cl*xi**2 + dl*xi**3 
mm 130 continue 

do 115 1=2,22 
do 115 j=2, 22 
x(i,J) = x(i,l) 

y(i.j) = (x(i,j)-x(l,j))/tan(del) 

115 continue 

do 120 i=2, 7 
1 do 120 j=2, 22 

z(i,j) = z(l,J) 


120 continue 

do 125 i=8, 22 
do 125 J=2, 22 

z(i,j) = z(i,l)*z(7,j)/z(7,2) 

125 continue 

if (k.eq.l) open (7,file='bluntlo.dba' ) 
if (k.eq.2) open (7, file=' interlo.dba' ) 
if (k.eq.3) open (7,file='sharplo.dba' ) 
if (k.eq.4) open (7, file=' extralo.dba' ) 
write (7,*) 22,22 

write (7,*) ((x(i, j)/rc,i=l,22), j=l,22), 

1 ((y(i.J)/rc,i=l,22),j=l,22), 

2 C(-z(i,J)/rc,i=l,22),j=l,22) 
close (7) 

900 continue 
end 
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program atinjoin 

dimension x(30, 10),z(30, 10), y( 30, 10) 

pi = 4.*atan( 1. ) 

del = 65 . /180.*pi 

re = 12.*tan(del) 

xte = . l*rc • 

data bl/. 510024/, cl/-. 19819008/, dl/. 025671542/ 
data as/. 3094608214/, bs/ .3327982281 /, cs/-. 04163681011/, 
1 ds/. 001507294514/ 
tmax = bl*xte + cl*xte**2 + dl*xte**3 
do 100 j=l,10 
x(l,J) = rc-10.022+.627 
y(i,J) = 0. 

z(l,j) ~ tmax 
do 120 i=2 , 15 
xi = (i-D/14. 

xi = .627 + (9.395-xte)*(exp(3.*xi)-l.)/(exp(3.)-l.) 
r = as*xi**.5 + bs*xi + cs*xi**2 + ds*xi**3 
t = tmax 

ang = (pi/2. -asin(t/r))/9.*( j-1) 
x(i, j) = rc-10. 022+xi 
y(i,j) = r*s in(ang) 
z(i,j) = r*cos(ang) 
if (j.eq.10) z(i,j) = tmax 
120 continue 

do 130 i=16, 30 
xi = (i-15)/15. 

xi = 10. 022-xte+xte*(exp(-4.*xi)-l. )/(exp(-4. )-l. ) 
r = as*xi**.5 + bs*xi + cs*xi**2 + ds*xi**3 
if (xi.gt.9.5) r=l. 65 
xt = 10.022-xi 

t = bl*xt + cl*xt**2 + dl*xt**3 
ang = (pi/2. -asin(t/r))/9.*( j-1) 
x(i,j) = rc - xt 
y(i» J) = r*sin(ang) 
z(i,j) = r*cos(ang) 

130 continue 
100 continue 

z(30, 10) = 0. 

open (7 , f ile=' jointup.dba' ) 
write (7,*) 30,10 

write (7,*) ((x(i, j)/rc,i=l,30), j=l,10), 

1 ((y(i,J)/rc,i=l, 30), j=l,10), 

2 ((z(i,j)/rc,i=l,30),j=l,10) 
close(7) 

end 


239 


program sting 

dimension x( 10) ,y(10, 10) ,z(10, 10) 

pi = 4.*atan(l. ) 

open (7,flle= , stlngup2.dba' ) 

del = 65./180.*pi 

rc = 12.*tan(del) 

x(l) = rc 

x(2) = rc + 4.5 

do 100 1=1,2 

do 100 j=l,10 

y(i, j) = 1. 65*sin(pi/18.*(j-l)) 
z(i,j) = 1.65*cos(pi/18.*(j-l)) 

100 continue 

do 200 1=3,6 
x(i) = x(i-l)+.5 

r = 52.57-(50.92**2-(.5*(i-2))**2)**.5 
do 200 j=l,10 

y(i,J) = r*sin(pi/18.*( j-1) ) 
z(i,j) = r*cos (pi/ 18 . *( j - 1) ) 

200 continue 

x(7) = rc + 17.59 
do 300 j=l,10 

y(7, J) = 2. 125*sin(pi/18.*(J-l)) 
z(7, j) = 2.125*cos(pi/18.*(J-l)) 

300 continue 

x(8) = rc + 19.5 
do 400 j=l,10 

y(8,J) = 2.875*sin(pi/18.*(j-l)) 
z(8, j) = 2.875*coa(pi/18.*( j-1)) 

400 continue 

x(9) = rc + 25. 7 
r = 2.875 + 6. 2*tan(2.*pi/180. ) 
do 450 j=l,10 

y(9, J) = r*sin(pi/18.*( j-1) ) 
z(9, J) = r*cos(pi/18.*( J -1) ) 

450 continue 

x(10) = 20.*rc 

r = 2.875 + (19.*rc - 19.5)*tan(2.*pl/180. ) 
do 500 j=l,10 

y(10, j) = r*s in(pl/ 18 . *( j - 1) ) 
z(10,J) = r*coB(pi/18.*(J-l)) 

500 continue 

do 600 1=1,10 
z(i, 10) = 0. 

600 continue 

write (7,*) 10,10 

write (7,*) ((x(i)/rc, 1=1, 10) , j=l, 10) , 


1 ((y(i» j)/rc, 1=1, 10) , j=l, 10) , ((z(i, J)/rc, 1=1, 10) , J-1, 10) 
close (7) 
end 
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Fig. 1 Wing and Sting Geometry (Right Half) 
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Fig. 2 Grid on Wing Surface 
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Fig. 3 Grid for Entire Computational Domain Around Wing Alone 
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Fig. 6 Blocking Strategy for Wing Model with Sting 
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Fig. 7 Seven-Block System for Entire Computational Domain 
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